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We present a comparative study of high carrier density transport in mono-, bi-, and trilayer 
graphene using electric- double- layer transistors to continuously tune the carrier density up to values 
exceeding 10 14 cm -2 . Whereas in monolayer the conductivity saturates, in bi- and trilayer filling of 
the higher energy bands is observed to cause a non-monotonic behavior of the conductivity, and a 
large increase in the quantum capacitance. These systematic trends not only show how the intrinsic 
high-density transport properties of graphene can be accessed by field-effect, but also demonstrate 
the robustness of ion-gated graphene, which is crucial for possible future applications. 

PACS numbers: Valid PACS appear here 



Transport through graphene is currently investigated 
in the low carrier density regime (n ~ 10 12 cm -2 ), where 
electrons behave as unusual chiral particles [TJ[2]. Despite 
exciting theoretical predictions (superconductivity [3H5]) 
and clear technological relevance (transparent electrodes 
[6], supercapacitors [7], and bio-sensors [8]), the high car- 
rier density regime has remained vastly unexplored due 
to the limited density range accessible in conventional 
transistors [TJ [2]. The recent development of so-called 
ionic-liquid gates, in which the coupling between gate 
electrode and transistor channel is realized through mov- 
ing ions that form an electric double layer (EDL) at the 
liquid/channel interface (Fig. la), is now changing the 
situation. With ionic-liquid gates, the gate voltage ap- 
plied up to several volts drops across the approximately 
1-nm thick EDL, leading to a very large geometrical ca- 
pacitance. As a result, the induced carrier density can 
easily exceed ti^d ~ 10 14 cm -2 , much larger than what 
is achieved in conventional solid-state field-effect transis- 
tors (FETs). Such a very strong field effect is valuable for 
technological applications (for instance, in organic FETs 
[9], where it enables low- voltage operation), and as a ver- 
satile and effective tool to tune electronic states in a rich 
variety of systems (by modulating metal insulator tran- 
sition [10], magnetoresistance [IT] , and by inducing su- 
perconductivity at the surface of insulators p~2j [13] ) . 

Recent work shows that ion gating can be used in com- 
bination with graphene. Experiments (e.g., Raman spec- 
troscopy [13] , quantum capacitance [15] , etc.) have fo- 
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FIG. 1: (a) A schematic cross-section of graphene EDLT, 
with the equivalent electrical circuit (b) A schematic repre- 
sentation of a device including the bias configuration used in 
the electrical measurements, (c) Optical microscope image of 
an actual device. One of two graphene device is immersed in 
the ionic liquid together with two Pt wires, acting as gate and 
quasi-reference electrodes. 



cused almost exclusively on specific properties of mono- 
layer at carrier density values up to ti^d ~ 5x 10 13 cm -2 , 
but no characteristic high carrier density features in the 
transport properties were identified. Here, as an efficient 
strategy to reveal these characteristic features, we per- 
form a comparative study of transport in ion-gated mono- 
, bi-, and trilayer graphene at high carrier density of 10 14 
cm -2 . The motivation for this strategy is twofold. First, 
when ri2D exceeds values of 10 13 cm -2 , differences be- 
tween monolayer and bi-/trilayer are expected, because 
in the latter systems the higher-energy split off bands 
start to be populated, which can provide an effective way 
to identify signatures of the intrinsic properties charac- 
teristic of the layers of different thickness. Second, in bi- 
and trilayer, the regime in which higher energy bands are 
populated has not yet been studied, and it is currently 
unknown how the opening of parallel transport channels 
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affects the conductivity of these systems. 

Mono-, bi-, and trilayer graphene devices were fabri- 
cated on Si02/Si substrates by exfoliating graphite (see 
Ref . [16] ) pQ . A small droplet of ionic liquid was applied 
onto the devices, covering the graphene layer, the gate, 
and a quasi-reference electrode, as shown in Fig. 1(b) 
and (c). The droplet can be easily removed and substi- 
tuted with a different ionic liquid, enabling the compari- 
son of EDLTs realized on the same graphene layer, with 
different ionic liquids (Fig. S2a in Ref. [16]): this is im- 
portant to check that the features observed in the experi- 
ments are not artifacts caused by the specific ionic liquid 
chosen. For all devices, the longitudinal sheet resistivity 
p xx , and Hall coefficient Ryl were simultaneously mea- 
sured at room temperature, in a Hall bar configuration, 
as a function of Vq — Kef ( Vg is the voltage applied on 
the Pt gate electrode and V re f is the voltage measured on 
the quasi-reference electrode as shown in Fig. 1(b). We 
found Kef ~ over the whole sweep range of Vq (Fig. 
S2b in Ref. [16]), which ensured almost all the applied Vq 
dropped at the liquid/graphene interface). The measure- 
ments were performed in a limited Vq range, to avoid the 
occurrence of chemical reactions between the ionic liquid 
and graphene, as it is necessary to obtain reproducible 
and reversible results (Fig. S3, S4 in Ref. [16 ). Despite 
this limitation, carrier density ~2x 10 14 cm -2 could be 
reached. 

Figure 2 shows the Vq dependence of the sheet conduc- 
tivity (j 2D = 1/Pxx for mono-, bi-, and trilayer graphene 
devices fabricated using ABIM-TFSI as ionic gate. For 
monolayer graphene (Fig. 2(a)), a linear increase of g^d 
is observed upon accumulating either electrons or holes, 
within AVq ~ =L1V from the charge neutrality point 
(CNP). For larger Vg, 0"2D exhibits a pronounced sat- 
uration [17]. The onset of a trend towards <J2D satura- 
tion is normally seen in conventional Si02-based mono- 
layer graphene FETs with sufficient high mobility [18]. 
Here, the use of EDLT makes the phenomenon unam- 
biguously clear, owing to the much larger carrier den- 
sity range spanned. The behavior of bi- and trilayer 
graphenes (see Fig. 2(b) and (c)) differs from that of 
monolayer. In particular, the linear increase of <T2D ap- 
pears within a narrower Vq range of A Vq ~ ±0.5 V near 
the CNP. More distinctly, outside this range the (J2D ex- 
hibits a non-monotonic behavior both for electrons and 
holes before keeping increasing further at higher Vq. 

Since charge accumulation in graphene EDLT devices 
is not simply described by a geometrical capacitance (see 
below), it is necessary to independently determine the 
sheet carrier density n^D as a function of Vg, in order to 
interpret the (J2D data. To this end we have simultane- 
ously measured the Hall resistance as a function of Vq. 
The blue lines in Fig. 2(d)-(f) display the resulting 
derived for devices fabricated on layers of different thick- 
ness. As expected, at the value of Vq corresponding to 
the (7 2D minimum n2B changes sign, confirming the shift 
of Fermi level E F across the CNP (Fig. S3 in Ref. [15]). 
Using ri2D determined from Hall measurements, we can 




FIG. 2: (a)-(c) <J2d( Vg) of mono-, bi-, and trilayer graphene 
devices realized using ABIM-TFSI ionic-liquid gates. The 
black dashed lines indicate the linearity of (J2d( Vg) near the 
CNP (vertical white lines), (d)-(f) ri2D as extracted from Hall 
measurements for graphene layers of different thickness, show- 
ing three perfectly overlapped data-sets measured at different 
magnetic field values (Fig. S3b in Ref. [E]). Plot (d)-(f) also 
show the capacitance C of the layers. In panel (e) and (f), 
the vertical arrows show coincidence in Vg of the anomalies 
in <72D and the increase in C. The insets of panel (a), (b), 
and (c) illustrate the band structure of mono-, bi-, and tri- 
layer graphene. The green shaded areas illustrate the range 
within which the Ef can be shifted. 
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FIG. 3: (a) <72D of bilayer graphene devices fabricated using 
three different ionic liquids as a function of ri2D (measured 
from the Hall effect). The inset illustrates the maximum value 
of U2D for electrons (white arrows) and for holes (black ar- 
rows) accessible with the different ionic liquids. The largest 
n 2 D reached in this study is 2 x 10 14 cm~ 2 using DEME-TFSI. 
(b) Dependence of ri2D on Vg (measured from the CNP) for 
mono-, bi-, and trilayer graphene devices, using DEME-TSFI 
as ionic liquid. The dashed black line is the sum of n2D of 
mono- and bilayer graphene, which compares well to the n2D 
measured in the trilayer. 



directly extract the mobility of the devices. We find max- 
imum values of 5.5, 3.5, and 9 x 10 3 cm 2 /Vs close to the 
neutrality point, for mono-, bi-, and trilayer graphene re- 
spectively; similar mobility values are also observed using 
other ionic liquid, for instance DEME-TFSI (Fig. S5 in 
Ref. [16]). 
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Having determined the n2r> of carriers, we extract the 
total capacitance of the devices from C = dn^D/dVQ. 
As shown in Fig. 2(d)-(f), when plotted as a function 
of Vq, a strong asymmetry between electrons and holes 
is clearly present, in sharp contrast to carrier accumula- 
tion using solid dielectrics. As electron-hole symmetry is 
known to approximately hold in graphene on the studied 
energy scale, we attribute the asymmetries observed to 
the properties of the EDLs. In large part, they originate 
from the different size of the positive and negative ions 
forming the ionic liquid, which are responsible for differ- 
ent thickness of the EDL for opposite polarities of the 
Vq. Indeed, the details of the asymmetry are different 
for different ionic liquids, as can be seen in Fig. 3(a) 
where the characteristics of a bilayer graphene device re- 
alized using ABIM-TFSI, DEME-TFSI, and EMIM-TFSI 
are shown (note also that when plotted as a function 
of ri2D, the asymmetry is considerably less pronounced). 
More importantly, however, the devices based on differ- 
ent liquids exhibit a fairly good agreement in the main 
features of the (J2D, including the absolute values, the 
non-monotonic behavior, and the position of the features 
as a function of n^D. This indicates that the features in 
the <72D are intrinsic to bilayer, which are not influenced 
by the specific ionic liquid used. Similar insensitivity to 
the ionic liquids was also observed in mono- and trilayer. 

It is apparent from Fig. 2 that C depends very strongly 
on Vq. This is because C, the total capacitance measured 
experimentally, is given by 1/(7 = 1/Cl + 1/Cq, where 
Cl is the geometrical capacitance of the EDL formed 
between Pt and graphene, and Cq is the so-called quan- 
tum capacitance associated to the finite density of states 
(DOS) of graphene. Owing to the large Cl of EDLs 
(several tens of /iF/cm 2 ) [15 , Cq dominates the total 
C, which is why the can be tuned by applying only 
small Vq (in conventional graphene FETs, Cq is nor- 
mally negligible, because the Cl is three orders of mag- 
nitude smaller than of EDLs typically, for a 300 nm Si02 
layer, 12 nF/cm -2 ). The dominant role of Cq natu- 
rally explains why the amount of carrier accumulated in 
graphene layers of different thickness (at the same Vq 
and using a same ionic liquid) is very different, as a di- 
rect consequence of the different DOS. 

The insets of Fig. 2(a)-(c) illustrate the main features 
of the band structures of graphene mono-, bi-, and tri- 
layer. Whereas in monolayer only two (valence and con- 
duction) linearly dispersing bands S2± touching at zero 
energy are present, bi- and trilayer graphenes have ad- 
ditional bands at higher energy, around 0.4 (2?2±) an d 
0.5-0.6 eV (T 2 ±) from the CNP, respectively (Sect. 3 
in Ref. [16]); filling of these bands occurs when the n2B 
exceeds ~ 2 and 7 x 10 13 cm -2 in these two cases (esti- 
mated from known band structure by neglecting changes 
induced by the perpendicular electric field generated 
by applied Vq). Indeed, the anomalies (i.e., the non- 
monotonic behavior) in the (J2d occur at density values, 
larger for trilayer than for bilayer, close to the ones esti- 
mated above. We therefore attribute the anomaly in <T2D 
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FIG. 4: (a) Measured data of <72d and C — dri2F> / cIVg as 
a function of ri2D for mono-, bi-, and trilayer graphene. The 
red and cyan pair of arrows point to the <72D anomaly (non- 
monotonic dependence of <J2d on t^d)- The dashed lines indi- 
cate how the (j 2D anomalies in bi- and trilayer, corresponds to 
the 7i2D where the capacitance increases, (b) Results of the- 
oretical simulations for the same quantities plotted in panel 
(a). To analyze the data theoretically, we have chose for the 
geometrical capacitance a realistic value of Cl = 40/iF/cm 2 . 



to the presence of a new scattering channel -inter-band 
scattering- that opens when the higher bands are filled, 
and that, as it is known from conventional semiconduc- 
tors [19 , can strongly reduce the carrier mobility. This 
is also why no anomaly is seen in the monolayer, since in 
that case no higher energy band is present. 

The interpretation of the <J2D anomalies in terms of 
inter-band scattering is further supported by the ob- 
served behavior of the capacitance of bi- and trilayer 
graphene (Fig. 2(e) and (f)), which exhibits a sharp in- 
crease concomitant with the <J2D anomalies since when 
the Fermi level enters the B<i± or T^± bands the DOS in- 
creases in a step-like manner, and so does the Cq (in the 
experiments the step-like increase is broadened by disor- 
der). Finally Fig. 3(b) shows the relation between ri2r> 
and Vq for mono-, bi-, and trilayer devices, from which 
we find that n2r> for the trilayer closely matches the sum 
of the densities in the mono- and bilayer. Such a relation 
is expected to approximately hold from the theoretical 
band structure of these systems, since the two lowest en- 
ergy bands in trilayer correspond approximately to the 
linear band of monolayer and the lowest energy quadratic 
band of bilayer (Sect. 3 in Ref. [16]). 

To substantiate all the considerations just made, we 
have analyzed the high-density electronic properties of 
mono-, bi-, and trilayer graphene in terms of a simple, 
well-defined theoretical mode [17]. As our only goal at 
this stage is to show that all the main features observed 
in the <T2D and C measurements naturally arise from a 
formally correct theoretical framework, for simplicity we 
have included only the effect of weak short-range scat- 
tering [17] and adopted a self-consistent Born approxi- 
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mation scheme [20] (for a precise quantitative analysis, 
other scattering mechanisms should be included as well, 
e.g., to capture the effect of long-range Coulomb poten- 
tial [21], resonant scattering [22], as well as ripples [23 J. 
The experimental data for devices realized using ABIM- 
TFSI are compared with theory in Fig. 4. At a qualita- 
tive level, the theory reproduces the trends seen in both 
0"2d(^2d) and C(ti2d) curves, for the graphene layers of 
different thickness. This includes the trend towards (J2D 
saturation in the monolayer, and the anomalous features 
in the cf^d and C in bi- and trilayer. In particular, the 
non-monotonic behavior of the calculated cf^d in bi- and 
trilayer can be traced back to the presence of interband 
scattering, which agrees with our initial interpretation 
and confirms the relevance of this process at high n^D. 
Note, finally, that in the model we have not included the 
modification of the band structure due to the perpendic- 
ular electric field generated by Vq (i.e. the opening of 
a gap in bilayer [24], and the modification of the band 
overlap in trilayer [25]); inclusion of these effects may 
lead to a better quantitative agreement for the values of 
ri2D at which higher energy bands are populated in bi- 
and trilayer. 

From the results obtained, we conclude that ionic- 
liquid gating is an effective and reliable technique to 



accumulate very large amounts of carriers in graphene- 
based materials. The 7&2D that have been achieved in this 
study are not far from those for which superconductiv- 
ity is observed in graphite intercalated compounds (e.g., 
KCg, ~ 5 x 10 14 cm -2 ) and, indeed, superconductivity 
in graphene at high ri2r> has been predicted theoretically 
j3H5]. It is likely that these higher values can be 
reached by biasing Vq over a broader range at lower tem- 
perature p~3j[26]. For applications, the high n2D leading 
to very high (J2D values, and the good compatibility of 
ionic-liquid gates with graphene are important for the re- 
alization of transparent electrodes in flat panel displays 
[6] , supercapacitors [7] and biosensors [8] . At the current 
stage, it appears that all the ingredients necessary for 
a rapid progress in directions of fundamental and tech- 
nological interest that combine the unique properties of 
graphene and ionic-liquid gates are already available. 
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